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I^ ■ The nonperturbative approach to soft high-energy hadron-hadron scattering, based on 

\Q • the analytic continuation of Wilson-loop correlation functions from Euclidean to Minkow- 



skian theory, allows to investigate the asymptotic energy dependence of hadron-hadron 



^ I total cross sections in lattice QCD. In this paper we will show, using best fits of the lattice 



data with proper functional forms satisfying unitarity and other physical constraints, how 
indications emerge in favor of a universal asymptotic high-energy behavior of the kind 
B log s for hadronic total cross sections. 
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1. Introduction 

The problem of predicting total cross sections at high energy from first principles 
is one of the oldest open problems of hadronic physics, not yet satisfactorily solved in 
QCD. Present-day experimental observations (up to a center-of-mass total energy ^/s = 
7 TeV, reached at the LHC pp collider P) seem to support the following asymptotic 
high-energy behavior: cr^^t i-^) ~ -Blog s, with a universal (i.e., not depending on the 
particular hadrons involved) coefficient B ~ 0.3 mb [2]. This behavior is consistent with 
the well-known Froissart-Lukaszuk-Martin (FLM) theorem [3], according to which, for 
s — 7- oo, o"tot (■5) < ^log^ (~)) where m^ is the pion mass and Sq is an unspecified 
squared mass scale. (Let us observe that the experimental value of B is much smaller 
than [about 0.5%] the coefficient -\ appearing in the FLM bound.) As we believe QCD 
to be the fundamental theory of strong interactions, we also expect that it correctly 
predicts from first principles the behavior of hadronic total cross sections with energy. 
Anyway, in spite of all the efforts, a satisfactory solution to this problem is still lacking. 
Theoretical supports to the universality of the coefficient B were found in the model of 
the iteration of soft-Pomeron exchanges by eikonal unitarization [1] (recently revisited in 
the context of holographic QCD [5]), and also using arguments based on the so-called 
Color Glass Condensate of QCD [6], or simply modifying the original Heisenberg's model 
[7] in connection with the presence of glueballs [8] . 

This problem is part of the more general problem of high-energy elastic scattering at 
low transferred momentum, the so-called soft high-energy scattering. As soft high-energy 
processes possess two different energy scales, the total center-of-mass energy squared 
s and the transferred momentum squared t, smaller than the typical energy scale of 
strong interactions (|t| < 1 GeV^ ^ s), we cannot fully rely on perturbation theory. A 
genuine nonperturbative approach in the framework of QCD has been proposed in [H] and 
further developed in a number of papers (see, e.g., [10] for a review and a complete list 
of references): using a functional integral approach, high-energy hadron-hadron elastic 
scattering amplitudes are shown to be governed by the correlation function (CF) of certain 
Wilson loops defined in Minkowski space [HI [121 [ISl [El HHj- Moreover, as it has been 
shown in [T6l [T71 [18], such a CF can be reconstructed by analytic continuation from 
the CF of two Euclidean Wilson loops, that can be calculated using the nonperturbative 
methods of Euclidean Field Theory. 

The analytic-continuation relations have allowed the nonperturbative investigation of 



correlators (and the corresponding scattering amplitudes) using some analytical models, 
such as the Stochastic Vacuum Model (SVM) [19], the Instanton Liquid Model (ILM) 
[20], [21], the AdS/CFT correspondence [22] and, finally, they have also allowed a numerical 
study by Monte Carlo simulations in Lattice Gauge Theory (LGT) [231 [21] (see also Refs. 
[23] for a short review). Although the numerical results obtained on the lattice can 
be considered "exact" (since they are derived from first principles of QCD), it is not 
possible to relate them directly to physical quantities, since the analytic continuation of 
the correlator can be performed only if an analytical dependence on the variables is known, 
while lattice data can be obtained only for a discrete (finite) set of values. However, it 
is possible to test the goodness of the known analytical models (like the ones we have 
mentioned above) simply through a best fit to the lattice data. This analysis has already 
been done in Refs. [23| [2T] and it is briefiy recalled in Section 3. The result of this analysis 
is not, generally speaking, satisfactory: known analytical models lead to bad quality best 
fits and, moreover, none of them provides a physically acceptable total cross section. 

In this paper, after a brief survey (for the benefit of the reader) of the nonperturbative 
approach to soft high-energy scattering in the case of meson-meson elastic scattering (in 
Section 2), and of the numerical approach based on LGT, comparing the numerical results 
to the existing analytical models (in Section 3), we will concentrate on the search for a new 
parameterization of the (Euclidean) correlator that, in order: i) fits well the lattice data; 
a) satisfies (after analytic continuation) the unitarity condition; and, most importantly. 
Hi) leads to a rising behavior of total cross sections at high energy, in agreement with 
experimental data. In particular, one is interested in the dependence of the CF on the 
angle 6 between the loops, since it is related, after analytic continuation, to the energy 
dependence of the scattering amplitudes, and also in its dependence on the impact- 
parameter distance. In Section 4 we show that, making some reasonable assumptions 
about the angular dependence and the impact-parameter dependence of the various terms 
in the parameterization, our approach leads quite "naturally" to total cross sections rising 
asymptotically as 5 log s (that is what experimental data seem to suggest). Moreover, 
in our approach the coefficient B turns out to be universal, i.e, the same for all hadronic 
scattering processes (as it also seems to be suggested by experimental data), being related 
to the mass-scale /j, which sets the large impact-parameter behavior of the correlator. This 
is actually the main result of this paper. In Section 5 we draw our conclusions and discuss 
some prospects for the future. 



2. High— energy meson— meson elastic scattering amplitude and 
Wilson— loop correlation functions 

We sketch here the nonperturbative approach to soft high-energy scattering (see [23] 
for a more detailed presentation). The elastic scattering amplitude Ai(hh) of two hadrons, 
or more precisely mesons (taken for simplicity with the same mass m), in the soft high- 
energy regime can be reconstructed, after folding with two proper squared hadron wave 
functions \ipi\'^ and \ip2\'^, describing the two interacting hadrons, from the scattering 
amplitude Ai^dd) of two dipoles of fixed transverse sizes -Ri,2±, and fixed longitudinal- 
momentum fractions /1.2 of the two quarks in the two dipoles [11], [T2|, [131 O [15] : 

M^hh)is,t) = jd^Rl^jdh\MRlLJl)\'[d'R2L[df2\MR2L,f2)\' 

X M^dd){s,t;Ri^,fi,R2±,f2), (2.1) 

with: / d^R,^ /; df, \MRi±, fi)? = I d'R2i. lo df2 \MR2±, /2)P = 1. 
(For the treatment of baryons, a similar, but more involved, picture can be adopted, us- 
ing a genuine three-body configuration or, alternatively and even more simply, a quark- 
diquark configuration: we refer the interested reader to the above-mentioned original 
references [IH dH [131 [H [15] .) 

In turn, the dipole-dipole (dd) scattering amplitude is obtained from the (properly nor- 
malized) CF of two Wilson loops in the fundamental representation, defined in Minkowski 
spacetime, running along the paths made up of the quark and antiquark classical straight- 
line trajectories, and thus forming a hyperbolic angle x — log(s/m^) in the longitudinal 
plane (see Fig. 1). The paths are cut at proper times ±T as an infrared regularization, 
and closed by straight-line "links" in the transverse plane, in order to ensure gauge in- 
variance. Eventually, the limit T — t- 00 has to be taken. It has been shown in [T6l [T71 [T8] 
that the relevant Minkowskian CF Qm{x'i T] z±] 1, 2) {z± being the impact param,eter, i.e., 
the transverse separation between the two dipoles) can be reconstructed, by means of 
analytic continuation, from the Euclidean CF of two Euclidean Wilson loops, 



-ig (h Af,{x)dXf, 

JCl,2 



(2.2) 
where (. . .)e is the average in the sense of the Euclidean QCD functional integral, and the 
arguments "1[2]" in Qe (and Qm) stand for "-Ri[2]±, /i[2]"- The Euclidean Wilson loops 



Wi2 are calculated on the following quark [q] - antiquark [q] straight-line paths, 

Ci : xf\T) = ~z + K + ffRu C, : Xfir) = K + f}^R,, (2.3) 



with r G [—T,T], and closed by straight-line paths in the transverse plane at r = ±T. 
The four-vectors pi and p2 are chosen to be pi_2 = ""^(isin |, 0_|_, cos |), 6 being the 
angle formed by the two trajectories, i.e., pi ■ P2 = m^cos6'. Moreover, Ri = (0,^j_l,0), 
z = (0, z±, 0) and // = 1 — /i, /f = — /j. We define also the Euclidean and Minkowskian 
CFs with the infrared cutoff removed as 



CE{e;zr,l,2) = \imgE{e;T-zr,l,2), 

T— >-oo 

Cm{x;z±;1,2) = lim ^m(x;T;£±;1,2). (2.4) 



T-i>oo 



The dd scattering amplitude is then obtained from CE{d', ■ ■ •) [with 6 G (0, vr)] by means 
of analytic continuation as 

Midd){s,t; 1,2) = -% 2s J d^z^e"^^-^^CM{x ^ log(s/m'); il; 1,2) 

^ -. 2./.-.e--C.(. ^ -., . -Uog(V™^);£.; 1.2). (2.5) 

where s and t = —\q±\'^ {q± being the transferred momentum) are the usual Mandelstam 
variables (for a detailed discussion on the analytic continuation see [H], where we have 
shown, on nonperturbative grounds, that the required analyticity hypotheses are indeed 
satisfied). By virtue of the optical theorem and of Eqs. (12.11) and (12. 5p . the total cross 
section is then given by the expression 

aL?(s) ~ -lmM^hh){s,t = 0) 



= -2jd^R,^J dh \MRi±Ji)\' J d'R2± J dh \MR2±J2)\' 

X f dh^ReCEiO ^ -IX ^ -ihgis/m^y, zr, Ri±, fu R2±j2)- (2.6) 

If one chooses hadron wave functions invariant under rotations and under the exchange 
/j — 7- 1 — /j (see Refs. [HI |15] and also jTO], §8.6, and references therein), the CF Ce in 
Eqs. (12.11) and (12. 6p can be substituted (without changing the result) with the following 



averaged CF [Ri± = |-Rj_l|(cos0j, sin^j)]: 

C^'"'(6'; |il|; |i?i_L|,/i, |-R2±|,/2) = / -^ / -^ 

Xt |C£;(6';£j_;_Ri_l,/i,-R2±,/2) +C£;(6';zl;i?i^, 1 — fi,R2±, f2) 

+ CE{e- il; /2i^, A, /22±, 1 - /2) + C,,(e; f^; i?i^, 1 - A, i?2±, 1 - A) } • (2.7) 

We note here that, as a consequence of the (Euchdean) crossing-symmetry relations |25j . 
C£;(7r — 6^; zj_; 1, 2) = C£;(6'; zl; 1, 2) = C£;(6'; Zi; 1, 2), where the arguments "z" stand for 
"— i?j_L, 1 — /j" {i = 1,2), the function C|f^ is automatically crossing-symmetric, i.e., 
C^^'^(7r — ^; . . .) = C'^^{9; . . .) for fixed values of the other variables. (The exchange "1, 2" 
— )■ "1,2", or "1,2" — )■ "1,2", as well a.s 6 ^ ir — 6, corresponds to the exchange from a 
loop-loop correlator to a loop-antiloop correlator, where an antiloop is obtained from a 
given loop by exchanging the quark and the antiquark trajectories.) 

3. Wilson— loop correlation functions on the lattice and 
comparison with known analytical results 

The gauge-invariant Wilson-loop CF Ce is a natural candidate for a lattice compu- 
tation. In Refs. [231 IS] a Monte Carlo calculation of Ce for several values of the relative 
angle and different configurations in the transverse plane has been performed, using 30000 
quenched configurations generated with the SU{3) Wilson action at /3 = 6/g'^ = 6.0, cor- 
responding to a lattice spacing a ~ 0.1 fm, on a 16^ hypercubic lattice with periodic 
boundary conditions. The Wilson-loop CFs have been constructed using loops with 
transverse sizes |f^i_L| = |r2±| = 1 in lattice units {Ri± = afi±, zj_ = adj_) and seven 
different values of the relative angle 6, i.e., cot 6* = 0, ±1, ±2. Without loss of generality 
(see the Appendix of Ref. [2T]), the longitudinal-momentum fractions have been taken 
to be /i = /2 = ^: the loop configurations in the transverse plane that have been studied 
are d± \\ ri± \\ f2± {''zz£^) and d± ± ri± \\ r2± C^^yy")- Also the orientation-averaged 
quantity ("awe") defined in Eq. (12. 7p has been measured. Finally, the CFs have been cal- 
culated for the values d = \d±\ = 0, 1, 2 of the transverse distance between the centers of 
the loops: as expected (see the discussion in Section 4.1 below), the CFs vanish rapidly as 
d increases, thus making a "brute-force" Monte Carlo calculation very difficult at larger 
distances. 



As already pointed out in the Introduction, numerical simulations of LGT can provide 
the Euclidean CF only for a finite set of 6'-values, and so its analytic properties can- 
not be directly attained; nevertheless, they are first-principles calculations that give us 
(within the errors) the true QCD expectation for this quantity. Approximate analytical 
calculations of this same CF have then to be compared with the lattice data, in order 
to test the goodness of the approximations involved. This can be done either by direct 
comparison, when a numerical prediction is available, or by fitting the lattice data with 
the functional form provided by a given model. The Euclidean CFs we are interested in 
have been evaluated in the Stochastic Vacuum Model (SVM) [19], in perturbation theory 
(PT) [261 [Til [I9], in the Instanton Liquid Model (ILM) [201E1], and, using the AdS/CFT 
correspondence, for the N" = 4 SYM theory at large Nc, large 't Hooft coupling and large 
distances between the loops [22], obtaining, respectively: 
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(0) 
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(3.1) 

(3.2) 
(3.3) 

(3.4) 



where the coefficients Ki = Ki{z±; 1,2) are functions of z± and of the dipole variables 
Ri±, fi- The comparison of the lattice data with these analytical calculations is not, 
generally speaking, fully satisfactory. The values of the chi-squared per degree of freedom 
(Xd.o.f.) of ^^^ best fits, performed in Ref. [23] using the above-reported functions (13. ip - 
(13. 4p . are hsted in Table [1] (together with the values obtained from best fits with the 
parameterizations "Corr 1" , "Corr 2" and "Corr 3" , that we shall introduce and discuss in 
the next section). As one can see in Table (H largely improved best fits have been obtained 
by combining the ILM and perturbative expressions into the following expression: 



smt^ 



(3.5) 



As we have said in the Introduction, the main motivation in studying soft high-energy 
scattering is that it can lead to a resolution of the total cross section puzzle. From this 
point of view, the analytical models considered in this section are absolutely unsatisfac- 
tory, since they do not lead to rising, or, better, to Froissart-like total cross sections of 
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0.55 
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0.15 
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3.89 


0.17 


0.16 
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0.12 
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Table 1: Chi-squared per degree of freedom for a best fit with the indicated function: 
SVM [Eq. dnH)], PT [Eq. (Q], ILM [Eq. ([SJD], ILMp [Eq. (133])], Corr 1 [Eq. (|4:23|) 
Corr 2 [Eq. g^S])], Corr 3 [Eq. I K29^ ]. 



the form Blog s at high energy, as experimental data seem to suggest. In fact, the SVM, 
PT, ILM and ILMp parameterizations f l3.ip -f l3l3|) and (13. 5p lead to asymptotically con- 
stant total cross sectionci, as it can be seen by using Eq. (12. 6p . Concerning the AdS/CFT 
expression (13. 4p . obtained in A^ = 4 SYM, it has been shown in [27] that, by combining 
the knowledge of the various coefficient functions Ki in (13. 4 p at large \z±\ [22] with the 
unitarity constraint in the small-|£i| region, a non-trivial high-energy behavior for the 
dd total cross section can emerge (including a Pomeron-like behavior a ~ s^^'^: note that, 
since a Conformal Field Theory has no mass gap, there is no need for the Froissart hound 
to hold also in this case). 



4. How a Froissart— like total cross section can be obtained: a 
new analysis of Wilson— loop correlators from lattice QCD 

An ambitious question that one can ask at this point is if the lattice data are compatible 
with rising total cross sections. An answer can in principle be obtained by performing best 
fits to the lattice data with more general functions, leading to a non-trivial dependence on 
energy. This approach requires special care, because of the analytic continuation necessary 
to obtain the physical amplitude from the Euclidean CF: one has therefore to restrict the 
set of admissible fitting functions by imposing physical constraints, ffist of all unitarity. 

Introducing the ^^ hadron-hadron correlator" Cj^j- as the average of the '^ dipole-dipole 
correlator" Cm over the dipole variables, weighted with the proper squared hadronic wave 



"Actually the ILM parameterization p.3p leads to null total cross sections! 



functions, i.e., 

cg'^^x;!^"!!) = J d'R,^ J dMMRi±Ji)\' J d'R2± J df2\MR2±j2)\' 

X Cm{x;z±;Ri^JuR2±J2), (4.1) 

one immediately recognizes from Eqs. (12. ip and (12. 5p (see Refs. [13] and [15]) that 
^M (XJ l^±l) is nothing but the scattering amphtude A{s, \z±\) in impact-parameter space 
(i.e., the partial-wave scattering amplitude Si — 1 = rjie^'^^^ ~ 1)) which must satisfy the 
following well-known unitarity condition (see, for example, Refs. [28]): \A{s, \z±\) + l\ < 1. 
Therefore, this unitarity condition immediately translates to: 

\ctt\xAz^\) + l\<l- (4.2) 

Since the hadronic wave functions are normalized to 1, the unitarity condition (14. 2 p is 
obviously satisfied if the following sufficient (and therefore ^^ stronger" ) condition for the 
loop-loop correlator Cm(x; ^±5 1) 2) holds: 

\CMix;Zr,Rl±JuR2±j2) + l\ <1 Vfx, Rl±, R2±, /l, /2, (4.3) 

i.e., if the dipole-dipole correlator stays inside the Argand circle for all values of z±, Rij_, 
R2±, fii /2- It is also possible to find another {^^weaker") sufficient unitarity condition, in 
terms of the averaged correlator C^'^, which is the Minkowskian version of Cf^'^, defined in 
Eq. (12. 7p in the Section 2: in fact, as we have said there, we can substitute the correlator 
Cm with the averaged correlator C^}"^ (without changing the result), whenever it is averaged 
over the dipole variables Ri_i and fi (with z = 1, 2) with the proper (squared) hadronic 
wave functions |^ip and |V^2p, as, for example, in Eq. (I4.ip . Therefore, the unitarity 
condition (14.20 is also satisfied if the following sufficient condition for the averaged loop- 
loop correlator {stronger than (14. 2p . but weaker than (14. 3p ) holds: 

|cr(x;l^±|;|i?i±l,/i,l^2±|,/2) + i|<i v|ii|, |i?i^|, |i?2±|, /i, /2. (4.4) 



4.1. General considerations on the form of the correlator 

In this section, we are going to introduce, and partially justify, new parameterizations 
of the CF that, in order: i) fit well the data; ii) satisfy the unitarity condition after 
analytic continuation; and iii) lead to total cross sections rising as 5 log s in the high- 
energy limit (as experimental data seem to suggest). 

9 



We will show that the above conditions lead to rather strong constraints about the possible 
shape of the parameterization. First of all, we observe that the conditions ii) and iii) 
cannot be simultaneously satisfied by a correlator with the following factorized form: 
—ReCj^ (x; |^±|) = '?(x)^(I^-l|)) with ^(x) rising with x for x — )■ oo, since the unitarity 
condition (14. 2 p implies that < — ReC]^ (xi l-^ll) ^ 2 (actually, with this factorized 
form, one can only have a^^^ (x) — > constant for x —^ oo). This means that rising total 
cross sections can be obtained, without violating unitarity, only if the correlator is not 
factorizable. Let us now give a few general considerations about the form of the Euclidean 
loop-loop correlator. As a starting point, we shall assume that the Euclidean correlator 
can be written as: 

CEi9; zr, 1, 2) = exp [Ks(^; z^, 1, 2)] - 1, (4.5) 

where Ke{0;z±; 1,2) is a real function (since the correlator Ce itself is known to be a 
real function [23]). This assumption, i.e., essentially the fact that: C^; + 1 > 0, is indeed 
rather well justified for many reasons: first, in the large-iVc expansion, the correlator Ce 
is expected to be of order 0{1/N^) (see Eq. (3.4) in Ref. [23]), so that C^^ + 1 > is 
certainly satisfied for large A^^; moreover, all the known analytical models (SVM, ILM, 
AdS/CFT correspondence, perturbation theory,. . . ) actually satisfy it; and last (but not 
least!), the lattice data obtained in Refs. [231 [2T] confirm ita 
At this point, the Minkowskian correlator can be obtained after analytic continuation: 



CMix; z±; 1, 2) = exp [KMix; z±; 1, 2)] - 1, (4.6) 

with /('^(x; ^±; 1, 2) = Ke{0 — )■ —ix'iZj^\l,2). In the large-x limit, the Minkowskian 
correlator Cm is expected to obey the unitarity condition (14. 3p . which using the parame- 
terization (14. 6p . reduces to the following very simple relation: 

Rei^M(x;^±;l,2)<0 \/z±,Ri^,fi (z = l,2). (4.7) 

So, the parameterizations that we are going to consider have the general form 

Ce{9; zr, 1, 2) = exp [ 5]; Ki{zr 1, 2)FEii9)\ - 1, (4.8) 



* Actually, lattice data [23] seem to support an even stronger condition, i.e., the positivity of the 
correlator itself, Ce > 0, which, in terms of the parameterization (|4.5|) would mean Ke > 0. We are 
not aware of any rigorous proof of this condition, apart from a rough argument in terms of "minimal 
surfaces" . (Let us observe also that this stronger condition is not expected to hold in the Abelian case: 
in fact, the parameterization for the Euclidean correlator Ce in quenched QED, found in Ref. |17| . while 
still being of the form (|4.5p . has not the property Ke > 0.) 

10 



where the sum is over different terms with various functions FEi{9) and various "param- 
eters" Ki{z±\ 1,2). We now make an important consideration about the dependence of 
the parameters Ki on the impact parameter \z_[_\. For a confining theory like QCD the 
loop-loop correlator Ce is expected to decay exponentially at large \z±] as 

Ci^-ae-'^l^^^', (4.9) 

where /i is some mass-scale proportional to the mass of the lightest glueball {Mq — 1.5 
GeV) or maybe (as suggested, for example, by the SVM model: see below) to the inverse 
iZ-^mc of the so-called vacuum correlation length X^ac-, which has been measured with 
Monte Carlo simulations on the lattice in Refs. [29] (see also Ref. [30] for a review), both 
in the quenched {X^ac — 0.22 fm) and full QCD {X^ac — 0.30 fm). (For example, in the 
SVM model, see Eq. (13?T1) . one finds [19] that, for \z±\ -)■ oo, Ksym ~ e"!^^!/^""", so that 
^SjVM ^ ^Kgyy^cot'^ 9 ~ ae~^l^-Ll, with /i = j^.) Therefore, we should require the same 
large-|ij_| behavior (14.91) for the parameters Ki in Eq. (14. 8p . i.e., Ki ~ e"'^'^^'. Instead, 
for a non-confining theory, e.g., for a conformal field theory, also different behaviors of 
the parameters K^ for large 1^11 are possible, typically like powers of 1/|^±|. This is 
what happens for the parameters -^^1,2,3 for the parameterization (13.41) obtained from the 
AdS/CFT correspondence [22], [27]. However, it can be shown that, in the general case, an 
asymptotic large- |^x| behavior of the parameters Ki hke powers of 1/|^±| leads to non- 
universal high-energy total cross sections, and can reproduce a Froissart-like behavior, 
(Ttot ~ B log^ s, only with very ''ad hoc" dependencies of the parameters K^ on powers 
of 1/1 ill and of the functions FEiiO) on powers of 6. 

4.2. How a Proissart like total cross section can be obtained 

Let us now assume that the leading term for x ~^ +c>o (i.e., for s — )■ 00) of the 
Minkowskian dipole-dipole CF is of the form 

Cm(x;^"1;1,2) ~ exp(2/3/(x)e-'^l^^^l)-l, (4.10) 

X-i>+oo 

where /3 = /3(1,2) is a function of the dipole variables and /(x) is a positive and real 
function rising with x, i.e, /(x) — > +C)0 for x ~^ +00: for example we can have fix) = 
e"'^, (coshx)", x^c"^, .... It is then clear that, in order to satisfy the unitarity condition 
( 14. 3p . the imaginary part of f3 has to be positive, i.e, 

|Cm(x;'?±;1,2) + 1|<1 ^ Im/3>0. (4.11) 

11 



The precise z^ dependence of fl4.10p is, of course, expected to be valid only for large 
enough |^j_|. For simplicity, we shall first assume that it is valid V|ij_| > 0. 
Following Eq. (12.61) . and performing the change of variable y = /i|^j_|, we have that: 

a'tJ:^ ^^Rej d^Ri^ j dh\MRi±Ji)\' J d'R2± J df2\MR2±j2)\'l{x, P), 

(4.12) 
where the quantity J(x, /?) is defined as 

J(x,/3)=/ dyy[l~exp{t/3f{x)e~y)]. (4.13) 

Jo 

We can now expand the exponential in series, and exchange the order of integration and 
summation, obtaining: 

;(,.,) . - f mMl rayye-', . _ f M^. (4,14) 



n! /n "' — ' nm'" 

n=l n=l 

The expression (14.140 is valid for an arbitrary function f{x)- However, as we have said, 
we are interested only in its asymptotic form for large Xi iii the case of f{x) rising with 
X, i-e., /(x) -^ +00 for X -^ +00. Thus 3xo e IR^ s.t. f\x) > 0, Vx > Xo- So, for 
X ^ X05 we can define the variable rj = log/(x), and re-write the quantity /(x, /?), Eqs. 
(I4.13p - (l4.14p . as a function of r] and /3: 

/(x,/3) = J(r/,/3) = -f;^^. (4.15) 

n=l 

Now, deriving (I4.15P with respect to the variable 77, one finds 

l^^-fW^. (4.16) 

ori ^-^ n\n 

The sum of the above series is known (see, e.g, Ref. [21]) and given by 

J'(r/,/3)=Ei(-z/3e'')+log(-z/3e^)+7, -vr < arg(-z/3e'') < vr, (4.17) 

where 7 is the Euler-Mascheroni constant (7 ~ 0.57721 . . . ) and -E'i(-z) is the Schloemilch's 
exponential integral (see, e.g., Ref. [31]). Since -£^1(2;) ~ e"^ / z at large \z\, for Re2; > 0, 
and, moreover, Re(— z/3e'') > <S=^ Im/3 > is nothing but the unitarity condition (14.111) . 
the asymptotic form of (14.171) is readily obtained: 

J'(r/,/3) ~ r/ + log(-^/3)+7 + C(e-''). (4.18) 

12 



Now we can re-integrate fl4.18p in 77, finding 

2 
J{r],(3) ~ — + 7/ (log(-z/3)+ 7) + constant + 0(e-''), (4.19) 

and, making the substitution 77 = log/(x), 

^(^' ^) ^^^00 ^^^Y^ + log /(X) ( log(-^/3) + 7) + constant + O (-^) . (4.20) 

Let us observe the following important fact: the leading order in x (and r]) does not 
depend on f3. So, coming back to Eq. (I4.12p . the asymptotic behavior of the total cross 
section turns out to be: 

^L? r^^Jd'Ri^ fdMMRi±Ji)\' Jd'R2± fdMMR2±j2)\' 
'iogV(x) 



+ log/(x)(log|/3|+7)+... 



(4.21) 



where the whole dependence on the dipole variables {Ri±, fi, i = 1,2) is in the coefficient 
p. If one assumes a leading term of the type f{x) = e"^ or, even more generally, f{x) = 
■^Pq^x in the correlator, the resulting asymptotic behavior for the total cross section is 
(recalling that x — log(s/m^)): 

^ff~i?log^5, with: 5 = ^. (4.22) 

We want to emphasize the fact that the above result is universal, depending only on the 
mass scale /x, which, as we have said, sets the large-|£i| dependence of the correlator. In 
fact, the integration over the dipole variables does not affect the coefficient of the leading 
term, since the hadronic wave functions are normalized to 1. Let us also observe, that the 
universal coefficient B is not affected by the masses of the scattering particles. In fact, 
in the case of scattering of two different mesons with masses mi and 7/12 ; the rapidity is 
X ~ log( — - — ). Therefore, considering the case of two different mesons (or, in general, a 
change of the energy scale implicitly contained in Eq. (I4.22p ) does not affect the universal 
coefficient B of the leading term (but will in general affect sub-leading log and constant 
terms), since x ~ log( — ' — ) = log(~) + log( '^° ), with ^sq being an arbitrary energy 
scale. This is the main theoretical achievement in this work. 

This same relation can be derived also with less stringent conditions on the z_[_ de- 
pendence, assuming (as it must be!) the exponential-type dependence in Eqs. (14. 9 p and 
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fl4.10p only for l^*^! > zq, with zq much larger than l//i and the dipole sizes. Starting 
from the expression fl2.6p of the total cross section, re-written as 

'^tot = ~47rRe/g°° |^±|c?|^±|Qf ix'i \^±\)^ where C\^ (x; \z±\) has been defined in Eq. 
(14. ip . one can split the integration in the variable \z±\ in two parts, a ^HaiF contribu- 
tion {J^ d\z±\ . . .), which can be evaluated using the approximate expression f l4.10p for 
the loop-loop correlator, and a ''core" contribution (J^° d\z±\ . . . ), which, instead, can be 
bounded using the unitarity condition (14. 2 p [or ( ]4.3p -( l44l) ]. in the form: —2 < Re Cm < 0. 

In the next subsection we are going to show our new analysis of the lattice data. 

4.3. New parameterizations for the correlator 

In what follows we show some parameterizations that we have found that satisfy the 
criteria i)-iii) listed above, together with the corresponding estimate of the asymptotic 
term of the high-energy total cross sections. Despite the appearances, it has not been 
simple to find such parameterizations: what follows is a selection between more than 
70 different parameterizations that we have tried. For each proposed parameterization 
of the correlator, we are going to show the Xd.o.f. ^^ ^^^ corresponding best fit to the 
lattice data for each given transverse distance and each given configuration: the results 
are summarized in Table [H As we have already pointed out in Section 2, the averaged 
correlator C^"^ is somehow "closer" to the hadron-hadron scattering matrix Ai(hh) than 
the correlator at fixed transverse configuration (like ''zzz" or ''zyy"), since it is actually 
the result of an integration of the dipole-dipole correlator over the orientations of the 
dipoles. The analysis performed above in Sections 4.1 and 4.2 can be repeated for the 
averaged correlator C"'^'^ without altering any conclusions. For this reason we are going 
to focus our analysis on the averaged correlator C™^ only; however, the Xdo f '^^ ^^^ t)est 
fits to the ''zzz" and the ''zyy" data are also shown for comparison in Table [1] Since, as 
we have already said in Section 2, the averaged correlator C"^ is automatically crossing- 
symmetric, so are the parameterizations that we proposeJH 



^However, it was observed in Ref. j23j that a small (but nonzero!) Odderon (C-odd) contribution in 
dd scattering, which is related through the crossing-symmetry relations |25| to the antisymmetric part 
of Ce{0) with respect to = -I, is present in the lattice data corresponding to the "zzz" and "zyy" 
transverse configurations. As noticed there, a crossing-antisymmetric term in the exponent Ke of the 
Euclidean CF (j4.5|) proportional to cot 9 (as, for example, the one appearing in the SVM parameterization 
(|3.ip and also the one appearing in the AdS/CFT parameterization p. 41) ) is in general suitable for taking 
into account this Odderon contribution and fits quite well the lattice data for the antisymmetric part of 
the correlator. Let us note also that its analytic continuation, coiO — > icothx, is limited for x — ^ oOj ^-nd 
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With regard to the exphcit angular dependence of the possible terms FEi{0), let us make 
some preliminary considerations. It is known [231 121] that lattice data for the correlator 
blow up at 6* = 0°, 180° (as expected from the relation between the correlator and the 
dipole-dipole static potential [32]), and that they are clearly different from zero for 6 = 
90°. A simple term like 1/sin^ (which always comes out as a Jacobian in the integration 
over the longitudinal coordinates in the analytical models considered in the previous 
section) can account for such a behavior and, as noticed in Refs. [23l [21], it fits quite 
well the data around 6 = 90°. Therefore, we shall always include a l/sin6' term in our 
parameterizations for the exponent in Eq. (14. 8p . 

Concerning the analysis of the impact-parameter dependence, we want to stress the fact 
that it must be taken only as an estimate, since only a few (small) values of the impact 
parameter are available from the lattice data (|^±| = ad, with d = 0,1,2). 

4.3.1. Correlator 1 

Following a first possible strategy, we have tried to improve best fits achieved with 
the ILMp expression (13.51) : the idea is to combine known QCD results and variations 
thereof. As an example, one could consider exponentiating the two-gluon exchange and 
the one-instanton contribution (i.e., the ILMp expression), and supplementing it with 
a term which could yield a rising cross section, e.g., a term proportional to cos ^ cot ^, 
like the one present in the AdS/CFT parameterization (13.41) . We thus find the following 
parameterization: 



Ce{9) = exp 



Ki , .. ,2 



sin 6* 
whose Minkowskian counterpart is: 



+ K2 cot^ O + K3 cos 6 cot 9 



(4.23) 



Ca/(x) = exp 



f -^1 , r^ „„„i, _ „„4-i, , A r.- „„iu2 



i — h K3 cosh X coth X ) — K2 coth x 

\ sinh X 



1. (4.24) 



The unitarity condition (14. 7p is satisfied if K2 > 0: from Table |2] one sees that the 
parameter K2 obtained from a best fit satisfies this condition, within the errors. The 
best-fit functions are plotted in Fig. [2l Performing a best fit with an exponential function 
~ e~^'^^' over the three distances, one finds that the coefficient K3 of the leading term 
for X -> 00 has a mass-scale fi = 4.64(2.38) GeV, that, following the result (I4.22p of the 



so it is consistent with the Pomeranchuk theorem, at least for rising total cross sections. 
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Corr 1 


d = 


d= 1 


d = 2 


K2 

Xd.o.f. 


5.85(42) ■ 10-3 

9.60(98) ■ 10-2 

-7.8(1.3) ■ 10-2 

2.81 


3.07(37) ■ 10-3 

2.44(49) ■ 10-2 

-1.37(72) ■ 10-2 

1.25 


8.7(3.1) -10-4 

-5.3(84.5) ■ 10-5 

1.7(1.9) • 10-3 

0.05 


Corr 2 


d = 


d= 1 


d = 2 


K2 

•x/2 
Xd.o.f. 


6.03(42) ■ 10-3 

4.63(46) ■ 10-1 

-4.54(50) ■ 10-1 

0.55 


3.26(38) ■ 10-3 

1.33(25) ■ 10-1 

-1.26(28) ■ 10-1 

0.31 


8.7(3.2) ■ 10-4 

-1.2(54.2) ■ 10-4 

1.7(6.7) ■ 10-3 

0.05 


Corr 3 


d = 


d= 1 


d = 2 


K2 

2 

Xd.o.i. 


6.02(36) • 10-3 

1.29(5) ■ 10-1 

0.17 


3.46(29) ■ 10-3 

4.47(27) ■ 10-2 

0.11 


1.07(20) ■ 10-3 

2.11(73)- 10-3 

0.10 



Table 2: Parameters (with their errors) for the Correlators 1 [Eq. (I4.23P ]. 2 [Eq. (I4.26P 
and 3 [Eq. (I4.29p ]. obtained from best fits to the averaged lattice data, and the corre- 
sponding Xd.o.f.' fo^ ^h^ transverse distances d = 0,1,2. 



previous section (with n = 1, as implied by this parameterization), leads to the following 
asymptotic total cross section: 



cr, 



(hh) 
tot 



5 log s, with: 



5 = 0.113+°;3^4 ^^ 



(4.25) 



This is compatible, within the large errors, with the experimental result B^xp — 0.3 mb 
reported in the Introduction. 

4.3.2. Correlator 2 



Another possible strategy is suggested again by the AdS/CFT expression (13. 4^ : one 
can try to adapt to the case of QCD the analytical expressions obtained in related models, 
such as A/" = 4 SYM. Although, of course, Eq. ( 13. 4p is not expected to describe QCD, it is 
sensible to assume in this case a similar functional form (basically assuming the existence 
of the yet unknown gravity dual for QCD). Assuming moreover that the known power-law 
behavior of the i^j's (expected for a conformal theory) goes over into an exponentially 
damped one (expected for a confining theory), K^ ~ e"'^'^^', one obtains a Froissart-like 

Blog^s. In this spirit, the second parameterization that we 



total cross section a. 



(hh) 
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propose is: 



Ce{9) = exp 



sin 6 



+ K2 



6* cot 6* + K-i cos e cot e 



(4.26) 



This one contains, in addition to the usual AdS/CFT-hke terms l/sin6', cot^ and 
cos 6* cot 6*, also another term proportional to 6' cot 6*0 the coefficients of the terms cot^ 
and 6 cot 6 are constrained by requiring that Ce{0) is crossing symmetric. The analytic 
continuation of fl4.26p is 



Cm{x) = exp 



sinhx 



vr 



fi'a-cothx 



K'^ cosh X coth X ] — X-^2 coth x 



1. (4.27) 



The unitarity condition (14. 7p becomes, in this case, K2 > 0, which is satisfied by the 
best-fit parameter within the errors (see Table [2]). The best-fit functions are plotted 
in Fig. [3l After the best fit over the distances with an exponential function, one finds 
for the leading-term coefficient K3 a mass-scale /i = 3.79(1.46) GeV. Thus, by virtue of 
Eq. (I4.22P (with n = 1, as implied by this parameterization), this correlator leads to the 
following asymptotic total cross section: 



a. 



{hh) 
tot 



B log^ s. 



with: B = O.nOtnil] mb 



0.081 



(4.28) 



that is again compatible, within the large errors, with the experimental result. 



4.3.3. Correlator 3 

The last parameterization that we are going to propose is 



Ce{0) = exp 



Binfl+^^Hs 



n 



9] cos9 



-1. 



(4.29) 



The ffist term is the usual l/sin6', while the second one is less "familiar", in the sense 
that it is not present in the analytical models known in the literature: but is a fact that, 
using this parameterization, the best ffi is extremely good (see Table 1), even if it has 
only two parameters. The Minkowskian version of the correlator (I4.29P is 

r3 



Ca/(x) = exp 



'■[ . -^ + K2 cosh X (tTT^X 
Vsmhx ^A 



X 



+ i\:2COshx(y--vrx^) 



1. (4.30) 



■i- Although such a term could seem "strange" at first sight, not being a "simple" combination of sin0 
and cos6', it has been shown in the first Ref. [TB], through an explicit calculation up to the order 0{g'^) 
in perturbation theory, that a similar term actually shows up in the case of the CF of two Wilson lines. 
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The unitarity condition f l4.7p reduces (in the large-x hmit) to K2 > 0, that is fully 
satisfied by the best-fit parameter shown in Table [2l The best-fit functions are plotted 
in Fig. m As regards the total cross section, let us note that in this case the leading 
term (for x ^ +C)o) in the exponent in (14.301) is of the form x^ ^^'- so, we can find the 
asymptotic behavior of the total cross section simply taking f{x) = X^ ^^ (i-6-, with our 
notation, n = 1 and p = 3) in the expression (I4.2ip . that leads again to the leading 
behavior reported in Eq. ( ]4.22p (with n = 1). 

After an exponential best fit over the distances, one finds that the mass-scale of the 
leading-term coefficient K2 is yU = 3.18(98) GeV. Thus, the asymptotic total cross section, 
derived from this correlator, reads: 

a£t''^~51og2s, with: 5 = 0.245+°:^^3 mb. (4.31) 

The comparison with the experimental asymptotic coefficient is extremely good and seems 
better than the previous ones, even if the errors are always very large. 

5. Conclusions 

The nonperturbative approach to soft high-energy hadron-hadron (dipole-dipole) 
scattering, based on the analytic continuation of Wilson-loop CFs from Euclidean to 
Minkowskian theory, makes possible the investigation of the problem of the asymptotic 
energy dependence of hadron-hadron total cross sections from the point of view of lattice 
QCD, by means of Monte Carlo numerical simulations. 

In this paper we have performed a new analysis of the data for the Wilson-loop 
correlator, originally obtained in Refs. [231 Elj by Monte Carlo simulations in Lattice 
Gauge Theory, and, in particular. Section 4 has been focused on the search for a new 
parameterization of the (Euclidean) correlator that, in order: i) fits well the lattice data; 
ii) satisfies (after analytic continuation) the unitarity condition; and, most importantly, 
in) leads to a rising behavior of total cross sections at high energy, in agreement with 
experimental data. In particular, one is interested in the dependence of the correlation 
function on the angle 6 between the loops, since it is related, after analytic continuation, 
to the energy dependence of the scattering amplitudes, and also in its dependence on the 
impact-parameter distance. In Section 4 we have shown that, making some reasonable 
assumptions about the angular dependence and the impact-parameter dependence of the 
various terms in the parameterization, our approach leads quite "naturally" to total cross 
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Corr 1 
Corr 2 
Corr 3 



/i (GeV) 
4.64(2.38) 
3.79(1.46) 

3.18(98) 



A = i (fm) 



rU.U45 
-0.014 

n ns2+°-°32 

U.U0Z„QQ^4 
0.028 
-0.015 



0.042: 



0.062 



B 



"2F" 

.,2 



:mb) 

hO.364 



0.113_oo37 

"-'•-'- '"-""0.081 
0.263 
0.100 



0.245' 



Table 3: Comparison of the mass-scale /i, the ^^ decay length'^ A 
B = 27r/yU^ derived from our parameterizations. 



l/fj, and the coefficient 



sections rising asymptotically as B log^ s (that is what experimental data seem to suggest). 
Moreover, in our approach the coefficient B turns out to be universal, i.e, the same for 
all hadronic scattering processes (as it also seems to be suggested by experimental data), 
being related (see Eq. 04.221) ) to the mass-scale fi which sets the large impact-parameter 
exponential behavior of the correlator: this type of behavior is typical of a confining 
theory, like QCD, and /i is expected to be proportional to the lightest glueball mass Mq 
or to the inverse of the so-called ^^vacuum correlation length" X^ac- This is actually the 
main result of this paper. 

Concerning the comparison between the numerical data obtained from the best fits and 
the experimental value of B, the agreement is quite good since the values are compatible 
within the large errors. In Table [3] we report the mass-scale /j, (and the decay length X = 
l//x) derived from the parameterizations that we have considered in the previous section, 
together with the predicted universal coefficient B = 27r//i^. However, we want to remark 
the fact that the values that we have found have been obtained from a limited set of "short" 
(i.e., surely not asymptotic!) distances, and so they must be taken only as an estimate. 
Of course, when more lattice data (at larger distances) will be available, the relation 
B = 27r/yU^ (derived from Eq. 04.221) with ri = 1, as occurs in our parameterizations of the 
correlator) may be confirmed or not. Concerning the relation between /i and the inverse of 
the vacuum correlation length X^ac or the lightest glueball mass Mq, at present a rigorous 
analytical determination in QCD is lacking (apart from the result fi = 2/Xyac obtained 
in the SVM model) and would be surely an important and helpful result. Using for fi an 
estimate derived from the experimental value of B, i.e, /ie^p = ^/'^n/Bexp — 2.85 GeV (as 
we have said above, our numerical estimates for /i are compatible with fi^xp within the 
large errors, as shown in Table [3]), one finds that fiexp ~ (3 ^ 4)/A^,ac or fiexp ~ 2Mg- Of 
course, only further investigations (both numerical and analytical) can confirm (or not) 
these results. In this respect, we must also remark that the whole analytic derivation of 
the result in Eq. (I4.22p is, of course, intended to be performed in full QCD (including 
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dynamical quarks), while the subsequent numerical analysis has been performed using the 
lattice data for the Wilson-loop correlation function which are available at the moment, 
and which were obtained in quenched (i.e., pure-gauge) QCD. However, we expect that 
the mass-scale //, which enters Eq. fl4.22p . is essentially gluonic, being related in some way 
(as we have said) to the vacuum correlation length X^ac or to the lightest glueball mass 
Mg, and so it should not dramatically change when including dynamical quark effects 
(hopefully, in near-future full-QCD lattice computations of the Wilson-loop correlator). 
Pushing this "speculation" a little bit further, we indeed expect, just on the basis of 
the experience with the vacuum correlation length X^ac (which increases from 0.22 fm in 
quenched QCD up to about 0.30 fm in full QCD), that the inclusion of dynamical quark 
effects should improve the agreement between the theoretical determination of fi (i.e., of 
the decay length X = l//i, and of the parameter B = 27r//i^: see the results in Table 
[3]) and its experimental value fi^xp — 2.85 GeV (corresponding to Xexp — 0.07 fm and 
Bexp — 0.3 mb), since A is expected to increase a little bit, so that /i should decrease and 
the parameter B should increase, moving towards the experimental value. Of course, it 
would be desirable to have lattice results in full QCD . . . 

Finally, let us observe that the functional integral approach turns out to be fundamen- 
tal for achieving this result: in fact, the investigation of hadron-hadron elastic scattering 
in the soft regime is mainly founded on the elementary loop-loop CF, which is then folded 
with some proper wave functions for the specific hadrons involved in the scattering pro- 
cess. Therefore, it is "natural" to expect that a universal behavior of the hadronic total 
cross sections at high energy may be originated by the loop-loop correlator itself: and 
actually it has been so. Of course, strictly speaking, our approach, based on the loop- 
loop CF, and the corresponding conclusion about the universality of B, only applies to 
meson-meson scattering: in this sense, we can consider as a real prediction the fact that 
the value of B that we have found from our analysis is consistent (within the errors) with 
the experimental value B^xp, which has been found considering baryon-baryon (mainly, 
pp and pp) and meson-baryon scattering. However, as briefly recalled at the beginning of 
Section 2, also for the treatment of baryons a similar, but more involved, picture can be 
adopted, using a genuine three-body configuration or, alternatively and even more simply, 
a quark-diquark configuration [HI [121 [131 [HI [15] . In particular, adopting a quark-diquark 
configuration for baryons, we can directly extend our approach (based on the loop-loop 
correlator) and the corresponding results to include also the case of baryon-baryon and 
meson-baryon scattering. This is probably enough to yield the Pomeron (C-even) contri- 
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bution (going as B log^ s) to hadron-hadron scattering, but not enough to consider also 
possible Odderon (C-odd) contributions, which are sub-leading in the high-energy limit, 
due to the Pomeranchuk theorem. In fact, as already noticed above (in Section 4.3), these 
C-odd contributions are averaged to zero in dipole-dipole scattering, since in this case the 
relevant CF Cf"*^ is automatically crossing-symmetric, and so they are probably visible 
only adopting a genuine three-body configuration for baryons. 
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Figure captions 



Fig. 1 The space-time configuration of tfie two Wilson loops Wi and W2 entering the ex- 
pression for the dipole-dipole elastic scattering amplitude in the high-energy limit. 

Fig. 2 Comparison of lattice data for the averaged correlator to best fits with the param- 
eterization fl4.23p (Correlator 1). 

Fig. 3 Comparison of lattice data for the averaged correlator to best fits with the param- 
eterization fl4.26p (Correlator 2). 

Fig. 4 Comparison of lattice data for the averaged correlator to best fits with the param- 
eterization (I4.29P (Correlator 3). 
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Figure 1: The space-time configuration of tlie two Wilson loops Wi and W2 entering the 
expression for the dipole-dipole elastic scattering amplitude in the high-energy limit. 
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Figure 2: Comparison of lattice data for the averaged correlator to best fits with the 
parameterization f l4.23p (Correlator 1). 
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Figure 3: Comparison of lattice data for the averaged correlator to best fits with the 
parameterization f l4.26p (Correlator 2). 
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Figure 4: Comparison of lattice data for the averaged correlator to best fits witli the 
parameterization f l4.29p (Correlator 3). 



